Solution of elastic-plastic shallow shell problems by the boundary element method

The boundary integral equations for elasto-plastic problems of shallow shells are established by using the fundamental solutions of shallow shells derived previously. The strains and stress-resultants in the plastic region are used as unknown variables. The simultaneous nonlinear equations of these variables and unknown boundary values are established and solved by direct iteration method.     

Vorticity shedding from a lentil-shaped body at large incidence in uniform flow

The 2D flow about a lentil-shaped body at high incidence in uniform stream is numerically simulated for high Reynolds numbers. The separation points are fixed on the body edges and the vorticity released through a Kutta-like condition is spatially discretized by point vortices. Their Lagrangianinviscid dynamics is described by using both aconformal mapping and aboundary elements method. The present results are in good agreement with previous calculations but still not satisfactory when compared with the experiments.The aim of the present work is to investigate the effect of different Kutta conditions and the effect of different amalgamation procedures in the far wake on global quantities such as theStrouhal number or the total shed circulation. Other remarkable effects induced either by the action of viscosity or by the rising of secondary 3D motion, may not be analyzed by the present simplified model.

---

Sommario Si studia il flusso attorno ad una lente ad elevata incidenza nelle ipotesi di elevato numero di Reynolds e flusso bidimensionale. I punti di separazione sono fissati sulle estremità della lente e la produzione di vorticità da questi è realizzata con opportune condizioni di Kutta. La dinamica dei vortici è descritta con un modello Lagrangiano non viscoso. I risultati, ottenuti sia con l'uso di una trasformazione conforme che con un metodo agli elementi di contorno, sono in ottimo accordo con i risultati di precedenti simulazioni numeriche, ma non sono ancora soddisfacenti se confrontati con gli esperimenti.Si analizza in termini di grandezze globali del flusso, quali il numero di Strouhal e la circolazione totale rilasciata, l'effetto di due diverse condizioni di Kutta e dell'uso di procedure di compattazione della scia lontana. Si conclude che l'uso di diverse condizioni di Kutta produce risultati pressoché identici mentre l'utilizzo di procedure di compattazione della scia può alterare in modo significativo il numero di Strouhal del flusso. Rimangono da analizzare altre possibili cause del disaccordo con le sperimentazioni quali la diffusione viscosa e la probabile presenza di moti tridimensionali.

     

奇性校正特解场法计算任意点应力和位移

本文将边界应力积分方程中的强奇性边界积分化为Ｃａｕｃｈｙ主值积分，证明了该主值积分的存在性，并给出它在变量替换中的比尺附加项显式；运用刚体位移法和单位应力场法将强奇性主子块和应力奇性系数矩阵表为同行副子块的线性组合，再配用极坐标变换和适当的数值求积计算各副子块。本法无论对平面问题还是空间问题，是光滑边界点、侧棱点还是角点，有限域或无限域，受载还是变温，都一概适用，而且可以扩展用来计算近边界点，公式统一，程序通用，数值效果良好。     

数值研究激波与旋涡的相互作用

从非定常形式的Euler方程出发,数值模拟了运动激波与旋涡相互作用的非定常流动过程。为保证激波具有较高的分辨率,采用对称型TVD格式进行了数值计算。结果表明。这样可以有效地模拟流场中一些复杂的流动现象,如激波变形、激波分叉和三波点的形成,以及旋涡结构的变化过程等,并与已有的实验流动显示相符良好。同时,也是对TVD格式求解这类问题的一次初步尝试。     

容量矩阵法在涡方法中的应用

本文将容量矩阵技术与快速Fourier变换相结合求解非矩形域的Poisson方程,再与格子涡方法相结合,发展了一套求解复杂外形的快速涡方法。作为算例,计算了有厚度平板和楔形体的分离流动,得到了满意的结果。特别是在垂直平板绕流中,模拟出了由于剪切层不稳定性引起的小涡结构。     

MACRO AND MICROSCOPIC PARAMETER MAPPING RELATIONSHIP FOR RED SANDSTONE 1)

颗粒流离散元软件PFC2D在岩土工程中的应用十分广泛,存在的主要问题是如何标定其细观参数,目前大多使用"试凑法",此方法的缺点在于工作量大、效率低等,本文提出以反向传播算法(back propagation,BP)神经网络的方式代替此方法,利用PFC2D 内置的FISH 以及Python 语言对其进行二次开发,使之自动运行并获取40组宏细观参数样本。结果表明：BP神经网络可以快速准确地建立宏细观参数映射关系,误差均保持在0.01 以内,模拟得到的应力--应变曲线与室内试验曲线高度吻合,且无需大量的数据样本便可创建网络,效率较高;另外,经验证本文选用的平直节理模型,可以有效地解决平行粘结模型UCS/TS值偏小问题,确定平直节理模型可以更好地模拟岩石。     

Tracking and reconstruction methods for moving-interfaces

A simple, efficient and accurate high resolution method to tracking moving-interfaces——the characteristic integral-averaging finite volume method on unstructured meshes is proposed. And some numerical tests and evaluation of six main efficient methods for interface reconstruction are made. Through strict numerical simulation, their characters, advantages and shortcomings are compared, analyzed and commended in particular.     

Disappearing Wakes and Dispersion in Numerically Simulated Flows Through Tube Bundles

Flow through a staggered array or bundle of parallel rigid cylinders of diameter D is computed with the help of a three-dimensional direct numerical simulation (DNS) at various values of Reynolds number between 50 and 6000. Two different spacings L of the tubes, i.e. L/D= 2 and L/D= 3, have been considered. When Re 500 the flow is laminar. In that case the converging flow between a pair of adjacent cylinders brings the oppositely signed vorticity at the two edges of the wake closer together behind the upstream cylinder so that the vorticity decreases quickly due to cancellation by diffusion. At Re 6000, when the flow is highly turbulent, the wake vorticity disappears rather by turbulent diffusion. This disappearance of the wakes in the closely packed flows (i.e. L/D 2) causes the mean flow in a cell, which consists of the region around a single cylinder, to be effectively independent of that in other cells. Another consequence is that the mean velocity field can be very well approximated by potential flow except in a thin boundary layer along the cylinder and a short wake behind it. The results have been applied to the transport of scalars in closely packed arrays. As in other complex flows, the dispersion of the scalars is dominated by the divergence and convergence of the streamlines around the cylinder rather than by the wake turbulence. Approximate expressions are derived for this topologically influenced dispersion in terms of the geometry of the array. The fact when most of the flow in the array can be approximated by a potential flow, allows us to introduce a fast approximate calculation method to compute the dispersion.     

功能梯度矩形厚板的三维热弹性分析

直接从三维热弹性力学基本方程出发，通过引入两个位移函数和两个应力函数，导出了一个二阶的齐次状态方程和一个四阶的非齐次状态方程。分析中采用了层合近似模型，即将板划分成厚度足够小的若干薄层，从面可将每一层内的材料常数近似为常数。给出了任意厚度的四边简支横观各向同性功能梯度矩形板的热弹性分析，特别当板较薄时，与薄板理论进行了数值对比，发现两者结果吻合很好。最后研究了材料梯度指标对热弹性场的影响，结果显示梯度指标对热应力和位移都有着显著的作用：在不同的区间，梯度指标对它们有不同的影响；并且在同一区间，梯度指标对两者的影响程度也有所不同。